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Abstract 



A new density matrix and corresponding quantum kinetic equations are introduced for 
fermions undergoing coherent evolution either in time (coherent particle production) or in 



O space (quantum reflection). A central element in our derivation is finding new spectral 

i_i solutions for the 2-point Green's functions written in the Wigner representation, that are 

,__! carrying the information of the quantum coherence. Physically observable density matrix is 

^ then defined from the bare singular 2-point function by convoluting it with the extrenous 

,-H information about the state of the system. The formalism is shown to reproduce familiar 

^ results from the Dirac equation approach, like Klein problem and nonlocal reflection from a 

1^ mass wall. The notion of the particle number in the presence of quantum coherence is shown 

O to be particularily transparent in the current picture. We extend the formalism to the case 

00 

(^ of mixing fields and show how the usual flavour mixing and oscillation of neutrinos emerges 

j> again from a singular shell structure. Finally, we show how the formalism can be extended 

*k> to include decohering interactions. 
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1 Introduction 

Many problems in modern particle physics and cosmology require setting up transport equa- 
tions for quantum systems in out-of-equilibrium conditions, including electroweak baryogen- 
esis [1], leptogenesis and particle creation in the early universe [3|, just to mention a few. 
Moreover, one often should be able to describe the evolution of coherent quantum correla- 
tions in the presence of decohering effects of the surroundings. This is true in particular for 
the case of the electroweak baryogensis (EWBG), where the problem is to reliably compute 
the fermionic out-of-equilibrium distribution functions in the neigbourhood of an expanding 
CP-violating phase transition wall. The EWBG problem can be divided to two regimes de- 
pending on the width of the phase transition wall in comparison to the mean free path of the 
fermions interacting with the wall. The case of thick wall has been addressed earlier in the 
semiclassical WKB approach in [H [S] (for earlier work see also [B] ) and later in the context 
of CTP formalism in [71 E]. For a review see [9]. In the thick wall limit the notion of local- 
ized particle distributions can be maintained and one can reduce the full quantum transport 
equations to local Boltzmann equations involving CP-violating (and CP-even) force terms 
employing a well defined expansion in spatial gradients. In the thin wall limit the dominant 
source for the asymmetry comes from the quantum reflection processes which are inherently 
nonlocal and no consistent quantum field theoretical formalism exists for treating reflection 
phenomena together with decohering collisions. For early attempts to include collisions in 
a Dirac equation approach see [IHl [HI [12]. One of the goals of this paper is to derive from 
field theory a density matrix formalism that can be used to solve this problem consistently. 
However, our methods can equally well be used to describe for exsample the coherent parti- 
cle production or neutrino-oscillations in the early universe. Although our motivation comes 
mostly from cosmological applications and we will only consider fermionic fields, the formal- 
ism that we will develop is not restricted to solving only these problems. Instead, our generic 
approach to coherence within quantum field theory should be easily extended to the case of 
scalar fields and also to nonrelativistic problems. It should have a wide range of applications 
in generic problems where one is interested in quantitative description of quantum coherence 
in noisy backgrounds. 

In this paper we shall consider only the noninteracting problem, but the formalism we 
develop is easily extendable to the case with interactions. Our main result is finding a phase 
space description for the quantum coherence in terms of new singular solutions in close 
analogy with the usual on-shell particle distributions. The basic objects of our study are the 
2-point Wightmann functionq^ 



^Note that our function G^ does not contain an explicit minus sign often included to its definition in the 
literature. 



iGaf3,ij{u,v) = (V'a,i(M)V^/3j(t;)) =Tr{p^«,i(M)^/3,i(f)}, (1) 

which describe the most interesting properties of the out-of-equihbrium fermionic system. In 
section |2] we will first derive the standard form of these functions under the usual assumption 
of translational invariance both in space and time and in thermal equilibrium. We then 
generalize this derivation to the case where the translational invariance is lost either in time 
or in one of the spatial directions (denoted by z hereafter). These studies are most easily done 
in a mixed representation, where the functions (fTl) are Fourier transformed with respect to 
the relative coordinate u — v. In this representation the equation of motion for G^ is found to 
separate into two sets of equations we call kinetic, or evolution equations containing explicit 
space or time derivatives, and to algebraic (in the mean field limit) cons^ramt equations. How 
this division takes place depends on the special assumption on the spacetime symmetries. 

Our most important observation is that giving up the translational invariance allows new 
solutions in the dynamical phase space that carry information on the quantum coherence in 
the system. In the time-dependent, but spatially homogenous case considered in section |3] 
this new class of solution is found at shell ko = and in the planar symmetric static case, 
studied in section El at shell kz = 0. (For a stationary problem, the position of the latter 
shell is shifted to k^ = v^ko, where v^, is the velocity of the static frame.) These solutions 
are interpreted to describe the coherence between particles and antiparticles of opposite 
3-momentas and same helicities on mass shells fco = ±v ^^ + m'^ in the former case and 
between left and right moving states of same spin on shells k^ = -^^jlz^^—rr^ in the latter 
case. The new coherence solutions are present only in the dynamical functions G*^ and G^ . 
We will show that the spectral sum rule excludes these solutions from the spectral function 
A = i{G^ + G^)/2, so that the coherence shells are not a part of the kinematical phase 
space. This is as it should be, since an asymptotic state made out of pure coherence without 
the mixing mass-shell states does not make any physical sense. Moreover, the coherence 
solutions in G^'^ are shown to be inconsistent both with the full translational invariance 
and with thermal equilibrium. 

The phase space shell structure described above is entirely set by the algebraic constraint 
equations. At first sight this singular structure appears to render the kinetic equations to be 
of little use, but in the end the problem reveals an interesting connection to the measurement 
theory. In section [5] we show how physical, observable density matrix can be defined as a 
convolution of the singular phase space solution with a smooth phase space weight function 
that describes the existing external information on the system. This definition of physical 
density matrices is the second major result in this paper, along with our finding of the 
singular coherence shells in the phase space, and we will illustrate this principle with several 
examples. We show, for example, how a complete information of the momentum, the energy 
and spin of the state renders the quantum evolution to a trivial constant propagation of an 
eigenstate without coherence. Other nontrivial examples with quantum mixing include the 



usual Klein problem, reflection off a smooth phase transition wall and the definition and 
evolution of the particle number in a homogenous out-of-equlibrium system. In all these 
cases we are able to reproduce the known results in a way which underlines the appropriate 
choice of the weight function and the necessity of including the new coherence solutions. 

In section [6] we extend our formalism to the case of mixing fields and show how the usual 
notion of the flavour mixing arises in the present context. Introducing flavour mixing through 
a Hermitian NxN mass matrix leads to a very complicated shell structure in the phase space. 
Assuming a suffiently large ko and no mass-degeneracies, each mass shell solution separates 
into A^ separate diagonal and A^(A^ — 1) off-diagonal mass shells. Similarly, the number of 
free coherence functions at shell kz = (in the static planar symmetric case) multiplies to A^ 
and in addition A^(A^ — 1) new coherence shells appear near, but not exactly at kz = 0. The 
role of each shell in the mixing phenomenon is described qualitatively. We show in particular 
how the usual density matrix equation for (flavour) mixing neutrino system emerges through 
a use of a weight function that encodes enough information to tell the direction of motion 
of the state, but not enough to collapse it into a singular mass-eigenstate. In section [7| we 
outline how our formalism can be straightforwardly extended to the case with interations 
and finally, section |8] contains our conclusions and outlook. 

2 Propagator theory 

In this section we will first review the standard derivation of a free fermion propagator in the 
thermal field theory, concentrating on the role of the underlying assumptions of translational 
invariance. We then extend the analysis to the case where the mass of the field can depend on 
the space and time coordinates. The loss of translational invariance leads to a rich structure 
of solutions for the 2-point function, and in particular to an emergence of quantum coherence 



as will be seen in sections 3.2 and 4.3 To be specific, we suppose that the Lagrangian of the 
theory is given by 

Ctee = ii'P^ + i'L'mi)R + i'Rm*tjJL, (2) 

where the mass m = m{t,x.) can be complex. This convention follows from the electroweak 
baryogenesis application where effective masses arising from higgs mechanism are spatially 
and temporally varying near the first order phase transition fronts and the complex mass 
is needed for CP-violation. Temporally varying mass term can arise for example in the 
context of particle prodution in the early universe [3]. For such nontrivial mass functions 
the Lagrangian £free can be understood as an effective theory for a fermion in a temporally 
and spatially varying background field. 



2.1 Standard derivation of a thermal propagator 

Let us first consider the case where m is constant. In the standard approach to thermal field 
theory in real time formalism, one introduces a complex time argument defined on some 
complex time path, an example of which is shown in figure [T| and introduces a propagator 
with complex time ordering 

Gc{u,v) = 9c{uq -vo)G^{u,v) - 9c{uo -vo)G'^{u,v) , (3) 

which, in the absence of interactions, obeys the equation: 

(z p^ -m*PL - mPR)Gc{u, v) = 6c{uo - Vo)6''{u - v) . (4) 

Here Pl,r = |(1 T 7^) and 9c{uo — vq) and Sc{uo — vq) are the step functon and the Dirac 
delta function on the complex time argument. It follows in particular that the Wightmann 
functions 

iG<{u,v) = {ij{v)ilj{u)) 

iG>{u,v) = {iIj{u)iIj{v)) (5) 

obey the homogenous equation: 

(ip^-m*PL-mPR)G<'>{u,v) = 0. (6) 

Let us now solve this equation under the assumption of translational invariance both in 
space and in time. Translational invariance is obviously only consistent with a constant m, 
and moreover it implies: 

G<>{u,v) = G<^>{u-v). (7) 

If we now define 

Gc (t, x) = - (z p +mPL + m*Pn) Ac (t, x) (8) 

it is easy to see that the functions A^'^ obey the Klein-Gordon equation 

{d^ - V' + |mnA<'> = 0. (9) 

Making a Fourier transformation w.r.t. the spatial coordinate, one finds the solution 

iA<'> (t, k) = al'^e''^'' + at'^e-'^''' , (10) 



where Uk = \J kP' + \'m\'^. The four independent coefficient functions (of k^) a^'^ can be 
solved in terms of, say a^, using the equation (6) and the identity A^{t,k) = A^{—t,k). 
The most general solution consistent with translational invariance is then: 



Im t li 



>► 

Re t 



Figure 1: Standard complex closed time path (CTP) contour for Keldysh propagators. 
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where a^ is the only yet unspecified function. If one assumes further that the system is in 
thermal equlibrium, one imposes the Kubo-Martin-Schwinger (KMS) boundary-condition: 



A>{t)^A<(t + tP). 



Using Eqs. (11 12) one then finds: 



1 



1+ = ^eql^fcj 



eP^k + 1 



and 



A> =e 
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where 



ZUk 



dko 



:i2i 



(13) 



(14) 



Using Eqs. (13 14) in (pj) one can write down the Mills representation for the full thermal 
propagator: 

(15) 



zGcit, k) = y ^ A{ko, k) Kq(A;o) + 0c(t)] e'""' 
where we have defined the spectral function 

A = 7rsgn(A;o)(^ + m^j — i'y^mf)6{k'^ — \m\'^) . 



(16) 



From the expression (16) one immediately finds the standard equilibrium propagators: 



zG< = 2Ancq{ko) and iG^ = 2^ (1 - n^q{ko)) . 



(17) 



Note that —2iA = G^ + Of as expected by the equal time anticommutator relations (see 



•eq 



section (3.3) below). We stress again that two crucial assumptions were necessary in arriving 



to the equilibrium expressions (15 17): the translational invariance both in space and in time 



and the standard KMS relation (12). 



2.2 Free fields in varying backgrounds 

Let us now assume that the mass is some nontrivial function of space and time m = m{t, x). 
Equation ^ still holds, but we can no longer use the translational invariance to simplify 
the problem. Instead, we can separate the dependence on the internal and external degrees 
of freedom by defining the Wigner-transformation as a Fourier transformation of a 2-point 
function w.r.t. the internal co-ordinate r = u — v. 

G{k, x)= I d^r e'^-'Gix + r/2, x - r/2) , (18) 

where x = {u + v)/2 is the average co-ordinate, and k is the internal momentum variable 
conjugate to u — v. Transforming ^ in this way gives: 

(H^ P. -^0 - tm,^')G<^^{k, x) = 0, (19) 

where rho and rhs are operators related to the real and imaginary parts of the mass function: 

^o,5G<'>(fc,x) = mj^j{x)e-'2^^-^^G<'>{k,x) . (20) 

Note that the derivative d^ operates only to the left, acting on the mass functions mR,i and 



the derivative d^ acts only to the right, on function G^'^ . Thus Eq. (19) contains an infinite 
number of derivative operators acting on functions G^'^. 

The dynamical functions G^'^ and the spectral function A are not Hermitian. It is 
desirable to work with Hermitian functions however, and to this end we will definejj 

G<'>{u,v)=iG<'>{u,v)^'^. (21) 

which are easily seen to be Hermitian in the sense that: 

[G^'^iu^v)]^ = G<'>{v,u) 

[G'<'>(A;,a;)]^ = G<'>{k,x). (22) 

It will also be convenient to write the equations of motion in the Weyl basis where the 
gamma-matrices are given by the following direct product expressions: 

7° = pi (g) 1 , a = -p^ (g) a and 7^ = -p^ 1 . (23) 

Here both p* and a* are the usual Pauli matrices such that the p-matrices refer to the chiral- 



and (T-matrices to the spin-degrees of freedom. In this representation, multiplying (19) from 
both sides by 7°, we find the equation: 

{ko + '-dt + p^(^[a-{k- '-V)] - (/mo - P^m,) ® I2) G<'>{k, x) = . (24) 

^Note that refs. [TJ |5] use a slightly different definition for the Hermitian function: G -^ i'j^G. Both 
definitions are equally correct, although they obey slightly different equations of motion. Our present 
convention is convenient in that iG^^^ is more directly related to the usual density matrix in Dirac indices. 



This is the master equation for this paper. It is difficult to analyse it further in full generality. 



However, Eq. (24) can be simplified by introducing certain space-time symmetries. The 
spatially homogenous case may be of interest to some applications in the early universe, 
when the spatial gradients may be neglected, but an evolution of some background field or 
the expansion of the universe creates nontrivial time dependence. This is the relevant limit 
for example for particle production in the early universe [3], coherent baryogenesis [13], or for 
the description of the neutrino-mixing in the early universe p!^. The static (or stationary) 
case with a planar symmetry is relevant for the electroweak baryogenesis studies. We will 
consider these special cases in sections|3]and|4]below, starting from the spatially homogenous 
time-dependent case. 

Before going to the special cases let us make the following genaralization to our master 
equation. As is well known, one can always consider the fundamental chiral fermions mass- 
less and introduce the masses through interactions, technically corresponding to insertion 
of the singular self-energy corrections. (This is actually the best way to understand the 
emergence of the spatially dependent mass terms.) However, we can also introduce interac- 
tions with other types of classical background fields by including the appropriate singular 



self-interaction terms. In this way we can us to generalize our master equation (19) to 



(H^P. -S,i„g(x)e~t9^9f)G'<.>(x,fc) = 0. (25) 

where Ssing(x) represents all relevant singular self-interaction corrections. We shall need this 
more general form when we consider the case of a quantum reflection from a potential wall in 
section 5.2 For now however, we will mostly concentrate to the case where SsiQg(x) reduces 



to the complex mass instertion ( 19 ) 



3 Spatially homogeneous case 

Let us ffist consider the case of a spatially homogenous system, where the translational 
invariance in space is restored. As a result the 3-momentum is conserved and the spatial 



gradient terms vanish in Eq. (24), giving rise to a much simpler equation for G^ 



(ko+'-dt~a-k- ^''rho - tiV^^) G<{k, t) = , (26) 



with 



^iar^E 



mo,,G<{k,t) = mR,i(x)e"^"* %G'<(A;,t) . (27) 



Of course G^ obeys an identical equation. Homogeneity also implies that helicity is a good 
quantum number. This can be seen from the fact that the helicity operator h = k ■ S = 



k ■ 7 77 , where k = k/\k\, commutes with the differential operator of Eq. (26). This 



fact is particularily transparent in the Weyl basis (23), where the helicity operator is just 



h = 1® k ■ a. As a result one can introduce a block- diagonal decomposition for G^ in the 
helicity basis: 



^h 



9h 



'l + hk- a) 



(28) 



where g^ are unknown 2x2 matrices in chiral indices. When the decomposition (28) is 
inserted into Eq. (26) it breaks into two independent equations (for h = ±1) for the g^- 
matrices: 

(ko + '-dt + /i|fc|p' - mop' + rh^p')gt = 0- (29) 

This is as far as one can simplify the equation by using the homogeneity. It is still useful 
to rewrite Eq. (29) in a form that separates the explicit dependence on dtg^. Taking the 
Hermitian and anti-Hermitian parts from (29) one finds two independent equations: 

2k,g< = Hg<+g<H^ 



idtgh 



Hg^ 



9hH^ 



(30) 
(31) 



where 



H = -h\k\p^ + rfiop^ - rh^p^ . (32) 

Note that the operators itiq 5, and therefore also the operator H are in general not Hermitian. 



Equations (30 31) receive a particularity nice interpretation in the mean field limit (where 
H becomes Hermitian), as will be described below. 

An alternative form, which we will find useful in our analysis, can be found by introducing 
the Bloch-representation for Hermitian g^: 



9t = 2 (^0 + i(^) 



(33) 



In this formulation g"^ is represented by a real 4-component vector instead of a Hermitian 
2x2 matrix. Using the Bloch-representation Eq. (29) can be written as: 





(^o+2^t)^o+ WAdl- ^0^1+ rn^di 



{ko + ^dt)g', + h\k\g, 



irhog!^ 



irh^Qi 



{ko + -dt)g^,-ih\k\gl^ 



rhogo + irhr^g^ 



{ko + -dt)g2 + ih\k\g'l + imog^ + ^5^0 





0. 



(34) 



These equations can again be separated into two independent equations by taking the real 
and imaginary parts. Formally these equations can be written as 

idtgl - "' ""^ 







(35) 
(36) 



where the matrix operators A^p and Ba^ are easily read off from Eq. (34). Equations (35) 
and (36), respectively, carry the same information as do the Equations (30) and (31). 



3.1 Mean field limit 

Transforming the original equation ^ for the two point function G^{u,v) to the mixed 



representation given by ( 18 ) resuhed in equations that contain arbitrary orders of derivative 



operators. As such it is difficuh to obtain any full solutions even after using the spatial 
homogeneity. To proceed further, we now consider the case where the gradient expansions 
are truncated to the zeroth order, i.e. the mean field limit. In this case the mass terms are 



no longer operators and the anti-Hermitian equation (30) becomes 



dtgh =-i[H.9h. 



while the Hermitian one (31) reduces to 

2kog^ = {H,g<}, 



(37) 



(38) 



where the mean field limit of the operator if. 



H 



-h\k\ 
m* 



m 

h\k\ 



(39) 



is immediately identified as the local Hamiltonian of the system. Anti-Hermitian component 
equations are often called kinetic equations (KE), whereas the Hermitian ones, which in 
the mean field limit are algebraic, are called '^constraint equations" (CE) [7]. Indeed, since 



equations (37 38) constitute 8 equations for 4 scalar quantities it is natural to interpret some 



of the equations as constraints on the phase space in which the dynamical solution is to be 
found [7j. 



The anti-Hermitian equation (37) looks very promising, since it clearly has just the stan- 
dard form of the equation of motion that one would derive for the density matrix ph = i^hi^ 
from the Dirac equation. However, it must be warned that interpreting it as a dynamical 
equation for g^, or even interpreting g^ as a density matrix, is not at all straightforward. 



Indeed, we will next find that as a result of constraint equations (38) the matrix g^ acquires 



a nontrivial singular structure, whereby the Eq. (37) is not even well defined as such 



3.2 Shell structure, homogenous case 



Let us first study the information contained in the constraint equations (38 ) in the mean field 



limit. The novel result of this section will be that the constraint equations allow, in addition 
to the usual free particle states, a class of apparently energy conservation breaking solutions. 
These solutions live on the shell /cq = and we interpret them as holding the information 
about the quantum coherence of mixing particle and antiparticle states (zitterbewegung). 
The complete shell structure imposed by the constraint equations is most easily seen by first 



rewriting them in the form (36), where Ba/d now is a constant matrix. We find 



kog^ + h\k\g^ - niRg^ + mig^ = 

ko93+h\k\g^ = 

kogi-niRg^ = 

ko92+mig^ = 0. (40) 

This set of equations has nontrivial solutions only when det^Bap) = 0. It is easy to see that 
in the homogenous case under investigation here, this condition becomes 

det(5„^) = kl{k^ - Imp) = (41) 

That is, in addition to the usual mass-shell solutions with k^ — |mp = 0, there are new 
solutions living on k^ = 0. Let us now find out the precise form of these solutions. 

3.2.1 ko ^ -solutions, mass-shell states 



From Eq. (40) it is easy to see that for k^ ^ the constraint equations (40) have the following 
solution: 

and 

(kl - \k\' - \mW, = Q. (43) 

This equation has the spectral solution: 

g^iko, \k\;t) = 27rf^^^i\klt) 5{k, - s,,u,) , (44) 

where Sk^ = sgn(/i;o) and the mass-shell energies are given by 



Uk = yk^ + |mp. (45) 



Using (42) and (44) we can write the full chiral mass-shell gf^-matrix as follows: 



,J„..(.„. |g|; t) ^ 2.|.„| /i;(|.1. (^ 1 ;J:l*l/*» ^ -^1^^^^ j .-(P - |„f ). (46) 

This solution has the expected form of a density matrix in the helicity eigenbasis. Exactly 
analogous solution exists for g^^_^- Invoking the usual Feynman-Stueckelberg interpretation 
the solutions with negative (s^,, = — 1) energies can be identified with antiparticles. Finally 
the two unknown functions f^{\k\,t) are generalized particle and antiparticle phase space 
densities. Note that in the limit m — i> we get the usual result that e.g. left chirality equals 
negative helicity for particles and positive helicity for antiparticles. 

10 



3.2.2 ko = -solutions, zitterbewegung 



The mass-shell solutions (46) were derived assuming that k^ ^ 0. However, setting ko = 



in the first place, but keeping \k\ ^ 0, we find out that equations (40) have a new class of 
solutions, which obey the relations 



9\ 



9'o 



^"^9^. 



u^I h 



(47) 



while the components §^2 are unconstrained. The corresponding spectral solution is 



^M-s(^0, |A?|;)^) = TT 




hmR/\k\ 



1 



1 —hmii/\k\ 

—hmj/\k\ —i 

i hini/\k\ 



5{k,) 



(4J 



where /{*(|fc|,t) and /^(|A;|,t) are two unknown real functions living on the shell k^ = 0, and 
so they cannot be directly associated with either particles or antiparticles. However, since 
one expects that a general time-dependent density matrix should contain information of the 
quantum coherence between particles and antiparticles, we make the following identification: 
the additional ko = ^-solutions (48) describe the quantum coherence (zitterbewegung) between 



particles and antiparticles with same helicity h and opposite momenta k. 

The most complete solution for a given momentum \k\ and helicity h can be written as 



9hiko,\k\;t) 



9h,m-siko,\k\;t) + g^Q_^{ko,\k\;t) 



(49) 



The full solution (49) contains /ottr independent spectral functions fa{\k\,t) living on three 



distinct shells. These shells are represented in the phase space diagram in figure [2j It is this 



singular structure which appears to render the evolution equation (37) to be of little use. 
This ambiguity is only lifted when one interprets g^ as a phase space density, and defines the 
true physical density matrix as a weighted integral over the singular g^. We shall postpone 
introducing these ideas more precisely until section [5j For now let us show that our new 
ko = 0-solutions are not present in the spectral function and in the usual thermal limits for 
G^'^ that follow by use of the KMS-relations. 



3.3 Spectral function 

The spectral function is defined in general as 



A{u,v) 



l{{Hu),i'{v)}) = liG>{u,v) + G<{u,v)) 



(50) 
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Figure 2: Dispersion relations in the case of a scalar mass function. The dark filled (red) 
blobs show the mass-shell contributions of a given \k\ to the matrix g'f^ and the light (yellow) 
blob shows the corresponding coherence solution from the new /cq = 0-shell. 

In the free field case A{u,v) obeys the same equation of motion as functions G^'^, and so 



the most general solution for A is of the form of Eq. ( 28 ) 



^7° = '^ah^-il + hk-a) 

h 



(51) 



where the chiral matrix a^ is identical to the most general solution (49) for g^, with four yet 



undefined spectral on-shell functions /^-^ for both helicities. However, in the homogenous 
case we are implicitly assuming that there are no nontrivial spatial correlations, so that we 
can perform the usual quantization by imposing the equal time anticommutation rules for 
the field operators: 

{i){t,u),i)\t,v)} = -i5^{u-v). (52) 

It is easy to see that these anticommutation relations imply that 



2^(t, u; t, {7)7° = 5^{u-v) 



(53) 



This condition is just the direct space version of the famous spectral sum-rule, whose mixed 
representation counterpart reads as: 



dfco 



TT 



^(fc,x)7° = l. 



(54) 



It turns out that the sum-rule (54) is enough to completely fix the values of the on-shell 
functions f^-^: 



f: 



hA 



and 



n 



hA 



0, 



(55) 



for both helicities. With these values the full solution for A becomes: 



A = 7rsgn(A;o)(^ + "^ij ~ i'y^mf)6{k'^ — \m\'^) . 



(56) 



12 



This is just the fainihar result for the spectral function in thermal quasiparticle limit found 
in section 2.1, see also for example 0. The spectral function is thus completely determined, 



and it does not contain any dynamics at all. Moreover, it does not have any contribution 
from the ko = 0-shell describing the coherence between particles and antiparticles. This is 
what one should expect: since coherence is a dynamic phenomenon, it should not show up 
in the measure of the one-particle phase space. This is precisely what we are seeing here. 
Moreover, it should be vanishing in the statistical equilibrium limit. 

3.4 Equilibrium limit for G^'^ 

The a priori independent distributions f^'^i, in functions G^ and G^ are constrained by the 
relation G^ + G^ = —2iA. Using Eq. (55) we then find: 



J!h< _j_ rh> 



fti< I fti. 
■1(1,2) "^ i(l 



'(1,2) 



h> 

(1,2) 



1 

0. 



(57) 



(Note that we drop the <, >-indices on /^-functions everywhere where there is no danger 



of confusion.) Relations (57) hold generally, as long as the spectral solutions are valid. 



However, if one further imposes thermal equilibrium by the Kubo-Martin-Schwinger (KMS) 
boundary condition 



Ge> (t) = Gf^it + 1(3) => G> (fco) = e^^°G< (fco) 
one can easily show that the quantum coherence functions must vanish: 



fh<,> _ r. 



(58) 



(59) 



and the mass-shell distributions become: 



/. 



h< 



£h> 

■'Sko 



neqiko) 

1 - neq(/Co) , 



(60) 



where neq(A;o) = l/(e^'^" + 1) is the usual Fermi-Dirac distribution. It is now easy to write 
down the equilibrium solutions for G^'^: 






/ . 5'eq,h 



[l + hk-a). 



(61) 



Using solutions (59 60) in (46), (48) and (49) equation (61) reduces to: 






27rsgn(A;o)(^ + mi?-z7^m7-)neq(/co)'^(^^ - \m\^) 
2']Tsgn{ko){ft + niR - i'y^mi) (1 - Ueqiko)) 6{P - \m\ 



(62) 
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which are just the standard equihbrium propagators found in section [2TT 

The crucial difference between our treatment and the standard derivation of the ther- 
mal propagator is how we treat the space-time symmetries. In the usual approach the 
coherence- solutions are excluded already before imposing the KMS-relations by the assump- 
tion of translational invariance. We can easily see how this limit arises in our approach. 



Eqs. (30 31) are already translationally invariant in space by homogeneity. Imposing also 



the time-translational invariance dtg^ = (and dtm = 0) turns the kinetic Eqs. (31) into 



four additional algebraic constraints. These constraints are consistent with the mass-shell 



solutions (46) as long as /^ are constants. Coherence solutions are inconsistent with them 
however, and hence the coherence is directly excluded by translational invariance. In our 
more complete treatment the functions fj} {\k\,t) can be time-dependent and the quantum 
coherence is maintained in the form of the dynamical functions fl[ 2){\k\,t). 

4 Planar symmetric case 

One often encounters situations where quantum states interact with classical backgrounds 
that can be approximated by planar configurations. Examples range from simple quantum 
reflection problems to particle interactions with an expanding phase transition wall during 
electroweak baryogenesis. Let us now assume that the system is symmetric along planes 
orthogonal to the z-axis. In this case the equation of motion (24) becomes 

(fco + ^dt - a' {K - ^5,) -a-k,- -f^nio - t-f'^^'m,^ G<{k;t,z) = 0. (63) 

Unlike in the homogenous case, helicity is not conserved here. However, one notices that 



apart from a,! ■ fcn-term the differential operator in Eq. (63) commutes with the spin in z- 
direction, which is described by the operator 5''^ = j^'j^'j^ = 1 (g) a^. One can try to get rid 
of the a||-terms by boosting to a frame where all reference to k^\ vanishes. Putting aside the 
transformation of the derivative-operators, the boost A,, should obviously be such that 

S{A,)^S~\A,) = kol''~k,^\ (64) 

where ko = sgn(ko)(/i;o — ^if)^- The explicit form of the boost matrix ^(Aii) is easy to work 
out: 

5(A||) = sgn(ko)^= ^ . (65) 

^J2~ko{ko + ~ko) 

The boost A,, obviously leaves the form of the derivative operator invariant: p = '-f^dt' — 
7ii ■ dx' — 'y^dz- However, as the boost mixes the time and space components, the planar 
symmetry argument 

d^^^G<ik;t,z) = (66) 
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does not hold anymore in the new coordinates. Instead, one can show that in the new 
coordinates 

a ■ ds> G'<{k; x') = 7,1 a ■ v^^ dtG'<{k; t, z), (67) 

where W|| = k^Jko and 7|| = 1/(1 — v^^Y^^. That is, the boost regenerates the noncomniuting 
a ■ ^1, -terms from the gradients and the boosted differential operator still fails to commute 
with S^ in general. An obvious exception to this rule is the static case where also dtG = 0. 



In the static case the boost (65) does remove all dependence on an from equation (63), and 
reduces the differential operator block-diagonal in the spin along the 2;-axis [7J. In what 
follows, we will restrict the analysis to the static, or more generally stationary cases. The 
latter can always be reduced to a static problem by a suitable Lorentz tranformation, as we 
shall see next. 

4.1 Stationary problems with planar symmetry 

In the application to the electroweak baryogenesis one can assume that the planar symmetric 
background fields have a stationary dependence on t and z. In particular for the mass function 
one can assume a form 

m{t,z) = miu{z — Vujt) , (68) 

where v^ is the velocity of the phase transition front in the plasma frame. This stationary 
form implies that the mass function is static in the wall frame, which is connected to the 
plasma frame by a Lorenz-transform 

tw = 1w{t - v^,z) ko^ = -fi,{ko - v^k^) 

Zw = lw{z - Vuit) ko^ = -fyj{k^ - Vyjkw) , (69) 

where 7^, = 1/(1 — f^)"*^/^. (That is: m(t,z) = mw{zw/'jw) = mizw)-) This boost is a 
constant in momentum variables, and hence leaves the mass-operators invariant. The spinor 
representation of the transform A^„ is 

5(A^) = -^(v/7- + l - ^/^^-l a""). (70) 

The boost S'(A^) obviously commutes with cJ||. In the wall frame the boosted function 

G<ik; zj = S{A^) G<ik; z - v^t) S{K) (71) 

obeys a static equation 

( fco„ - «=^ {K^ - ^9,J - « ■ fc|, - 7°mo - 27°7'^5 ) G<{k- z^)=Q. (72) 

where in particular the mass operators are static in the wall frame variables: 

"^0,5 = mR,i(2;^)e^^^"^^£» . (73) 
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Because the problem is static in the wall frame, Eq. (72), it can be boosted to a frame where 



all a ■ fc||-terms vanish. The explicit spinor representation ^(Aii^u) of the boost is given by 

(74) 



Eq. (65) with k^ -^ ko^. After this second boost one finds that the function 



^wwi^OiuJ kzn,'i ^w) = S{A\\^)G^{k; z^)S'(A||^) 



obeys the equation, 

ko^ -a^{kz 



-dzj - 7°^o - i7^7^m5 ) G<{kQ^,k-,^]z^) = 



(75) 



where ko^ = sgn{ko){kQ^ — k'^^Y^'^. In the doubly boosted frame, one can make use of the 
commutativity of the differential operator with S*^ and introduce the spin-decomposition 



analogous to Eq. ( 28 ) : 



l|Wv 



1 



9un.,s^2^^ + scTz), (76) 

where g^^ ^ are two (for s = ±1) unknown Hermitian 2x2 matrices in chiral indices. Inserting 



(76) into Eq. (^75^ one finds 



(ko^ + s{k,^ - -d,Jp^ - rhop^ + m^p'^jg^^^^ = . 



(77) 



This equation is very similar to Eq. (29). If one identifies sk^^ with h\k\ (and ko^ with ko), 
the only difference is replacing t^roperator by an operator —sp^dz- This small change leads 
to profoundly different solutions however. 

4.2 Division to constraint and evolution equations 



We now proceed to analyse (77) in a same manner as we analysed Eq. (29) in section 3 



However, a direct division of Eq. (77) into Hermitian and anti-Hermitian equations does not 



lead to the desired separation to kinetic and constraint equations Q Instead, one has to first 
multiply the equation from left by p^ and only then perform the division. In this way we 
find the equations 



2skzgf 
isd^gf 



Pgf + gfP^ 



where 



P = kop^ + iirhop"^ + m^p^] 



(78) 
(79) 

(80) 



''Note that in the earlier work related to electroweak baryogenesis |7] the conceptually wrong division to 
constraint and kinetic equations - following Hermiticity properties - was used. While certainly wrong for 
the discussion of the quantum reflection, we have checked that this choice in the end does not affect the 
semiclassical limit discussed in these papers. 
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Here we have dropped all indices referring to wall frame or to the frame with zero parallel 
momentum. One can either assume that boosts have been done, or that we consider the 
case with either t;^ = or A;|| = 0, or both. Obviously, the operator P is a generalization 
of a local momentum operator in the same manner as H generalized the local Hamiltonian. 
Alternatively, introducing again the Bloch representation 

9f^\{9o + 9tf^) , (81) 



we can decompose equation (77) into components as follows: 



1 
s{k^ - -d^)gl + kogl - irhog^ - irh^gl = 

s{.kz - -dz)g{ - ikog'^ + mofl's - i^59o = 

s{kz - ^d,)g^ + ikogl - imog'^ - m^g'^ = . (82) 

Again, these equations can be divided to real and imaginary parts, resulting in vector equa- 
tions isdzg^ = ^aiiQp and = B^^g^^, where A^^ and -B^^ are operators that can easily be 



read off from Eq. (82) and which become simple constant matrices in the mean field limit. 

Unlike H the operator P is not Hermitian even in the mean field limit. It is easy to see 
that it nevertheless takes the familiar form of the local momentum operator 

^-(:.::). («3) 

while the constraint and evolution equations become: 

-2skzgf = Pgf+gfP^ (84) 

isdzgf = Pgf-gfP^. (85) 

One again recognizes that the evolution equation has the standard form of the equation for 
a density matrix p = ipip'^ , which can be simply derived using the static Dirac equation [15j . 

4.3 Shell structure, planar case 

The analysis of the shell structure proceeds very similarly to the homogenous case. The 
new result will be that the constraint equations allow a class of apparently momentum 
conservation breaking solutions in addition to the usual free particle states. These solutions 
live on the shell A;^ = and turn out to hold the information about the quantum coherence 
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between mixing opposite momentum states with same spin. The shell structure becomes 
evident by solving the determinant condition for the mean field limit constraint equations: 



+ sk^gl - magi + mig^ = 

Ms + sMo = 
sk^gl + irLRgl = 
shg'^-migl = 0. (86) 

It is easy to see that the determinant of this set of equations is just 

det{B^^,) = kliP - \m\') . (87) 

So, by setting the determinant to zero we find again similar on-shell solutions as we did in 
the homogenous case, but instead of fco = solution, we now find a double root at kz = 0. 
Let us now find the explicit gf matrices corresponding to these solutions. 

4.3.1 kz ^ -solutions; free particle-shells 



It is again easy to show that for k^ j^ equations (86) have the particlular solution 



<! skz , , rriR „ mj , ,„„, 

93 = -T^9o, 9{ = ^gl g2 = -^g'o. (88) 

Kq Kq Kq 

and 

{kl-kl-\m\')gl = Q. (89) 



Equation (89) has the spectral solution 



\k I 
^^(/co, kz] z) = 271 fl^^ {ko, z) [-^ 5{kz - Sk^km) , (90) 

where s^^ = sgn(A;2). These solutions thus live on a well defined energy-momentum shell 
corresponding to the usual dispersion relation 



km = ^kl - |m|2. (91) 

Note that we are taking energy as the free variable, whereas kz is defined by the on-shell 
condition. This is the appropriate choice for a problem with spatial gradients, where the 



momentum need not be conserved. Using (88) and (90) we can write corresponding full 
chiral gf^-matrix as follows: 

S,<„_.(«=„. fc; z) = 2,|*„| /i . (*=„. .) ( 1 -ffl^" ^ "/J-)^. ) 5(e - \mn (92) 

This solution again has the expected form of an mass-shell state of a definite spin in z- 
direction, and again, an analogous solution exists for g^^_s- The unknown functions f±{ko, z) 
will be identified as generalized phase space densities for right and left moving states. 



4.3.2 kz = 0-solutions and quantum coherence 



The free particle solutions (92) were derived assuming that kz 7^ 0. Giving up this restriction, 



we find that equations (86) allow a new class of solutions living on shell kz = 0. Still keeping 



/co 7^ 0, we find the solution: 

rriR ^ mj ^ 

93 = 0, (93) 

while the components gl 2 are unconstrained. The corresponding spectral solution is 

gf^o~s{ko,kz;z) = TT f^{ko,z)i^^^ ° mn/ko J 

where /f and /| are unknown functions that only depend on the energy and the position. 
To see what physics these new solutions describe, note first that each state of a definite 
energy ko and spin s can correspond to two different states with opposite helicities and 
momenta. Secondly, a density matrix should carry information about the quantum coherence 
between degenerate states, that may be present when the defining quantum numbers (here 
the momentum) are sufficiently poorly known. No combination of free particle solutions 



(92) can carry such information however. We thus make the following interpretation: the 



additional kz=0-shell solutions (95) describe the quantum coherence of states of equal spin 
travelling in opposite directions. This is also very natural from the momentum conservation 
point of view: while the mixing mass shell components have large and opposite momenta 
kz = ±a/A;q — m"^, their coherent mixture has the momentum expectation value of kz = 0. 



The most complete solution satisfying the constraint equations (86) for a given spin s and 
energy /cq 7^ is 

gf{ko, kz] z) = ^<„,_s(fco, kz] z) + ^<o„s(^o, kz\ z) . (95) 

The practical use of this solution is again restricted by the fact that the matrices describing 
the physical mass-shell solutions and their quantum coherence are proportional to distribu- 
tions defined on different momentum shells in the phase space. These dispersion relations 
with particluar solutions for a given /cq are shown in figure [3j We will discuss the physical 
interpretation of these spectral solutions in section [5| Before entering this discussion we 
wish to end this section by showing how the coherence solutions change if one considers a 
stationary instead of a static problem: that is if one considers the problem in the plasma- 
rather than in the wall frame. 
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Figure 3: Dispersion relation in the case of a scalar mass term. The filled blobs show the 
on-shell contributions to the matrix gf. Peaked functions on each axis illustrate the accuracy 
of the exterior knowledge on fco and kz- In the case of kz the dashed line function illustrates 
an effectively accurate measurement of momentum direction and the solid line function the 
case of a complete ignorance on the direction. 

4.4 Finite v^; stationary case in wall frame 

Since the purpose of this section is merely to illustrate the effect of stationarity on the kz = 0- 
solution, we shall restrict ourselves to the 1+1-dimensional case here {i.e. we set ^n = 0). In 
that case the relevant equations of motion can be read from (63) with constant mass terms 
and replacement dt — > —Vwdz- Following our division to the constraint and kinetic equations 
in the static case we now find 



and 



fcofi-o + skzg's - mngl + rmg^ = 

kog^ + skzgo = 

-Vwdzg'l - 2.skzg2 + 2"^/5'3 = 

-Vwdzg2 + '^sk^gl + 2mRgl = , 

sdzg's = -v^dzgQ 

sdzgl = -^niRgl - 2mig{ - v^^dzgl 

sdzg{ = -2kog2 - 2mig'o 

sdzgi = 2kogl - 2mRgl . 



(96) 



(97) 



At first sight the situation appears problematic, since we now have only two purely algebraic 



constraint equations left in (96). However, since we have four independent linear equations 



for the two derivatives 5^5'i2; we can use the (last two) kinetic equations (97) to eliminate 



these derivatives from constraints. When this is done, the last two equations in (96) become: 

s{kz - v^ko)gl + mR{gl + sv^g^) = 
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sik^-Vu,ko)g2-mi{g'3 + sv^gQ) = 0. (98) 

We thus are back to four independent algebraic constraint equations. Using the second 



constraint in (96), which still gives (when fco ¥" ^) 93 = ~^9o^ the equations (98) can be 



rewritten as 

siK-vMigl + '^gi) = 0. (99) 

The presence of the new coherence shell solution is in fact more apparent here than it was 



in the static case: from (99) one immediately sees that if k^ — Vyjk^ ^ then gl = '^g^ 
and f^l = — y^gfg, which immediately leads to the usual mass-shell solutions. However, in 
addition the combinations gf = gf — j^g^ and ^2 = 5'2 + ^9o have spectral solutions 
^12 0^ fi2^ikz — Vwko). While the new shell structure is most obviously visible in variables 
gl 2 it is a simple matter of linear mapping to replace the free on-shell functions fl 2 with 
those associated with the distributions gf 2- The coherence shell has thus moved to k^ = w^ko 
as it should by the Lorentz transformation between the static and stationary frames. Setting 
kz = VyjkQ in the first two constraint equations then gives the g^ 3 in terms of the free variables 
gl 2 on the coherence shell: 



.s _ „., „ _, s_ 2,'^R s , ^I 



g'3 = -sv^go and ^o = -7^(^^i + ^^2) ■ (100) 

Kq Kq 

Finally, a direct calculation shows that the physical mass shells remain to be given by the 
Lorentz- invariant relation /cq — fc^ — |mp = 0. 

5 Dynamical equations and connection to the measure- 
ment theory 

In previous sections we have discovered and interpreted the complete shell structure of the 
free fermion propagator when a complete translational invariance is lifted either in the tem- 
poral, or in one of the spatial directions. We found in particular that in the mean field limit 
the propagator matrix has spectral solutions including the usual mass shell, but moreover 
also new solutions on shells where k^ = (homogenous) or fc^ = (static, planar symmetry). 
We interpreted these shells as carrying information about the quantum coherence between 
particles and antiparticles of equal helicity and opposite momenta in the homogenous case, 
and between left and right moving states of equal spin under reflection from a wall in the 
planar symmetric case. However, these interpretations are still problematic in that it is not 
clear what we mean by coherence, since it is living on a different singular shell. Indeed, both 
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our final dynamical equations (37) and (85) are yet ill defined because of this issue. We now 
show how this situation is to be interpreted, and in process discover an intersting connection 
to the measurement theory. 

5.1 Prom g'^ to a weighted density matrix 

The key idea is that in reality we can never have a complete information on the variables 
describing a certain process. Hence the physically interesting object - whose evolution we 
can study - is not the singular density matrix gf^^ but some smeared-out objedrl whose 
definition involves the extrenous information about the parameters of the system into the 
theory. To see how this works consider first the planar symmetric case, where we assume 
that we have a fairly precise information of the momentum k^ as well as of the energy k^ and 
the spin s of the state at all z. This situation is illustrated in the figure [3] for the phase space 



variables. This information excludes the coherence solutions and the full density matrix (95) 



is reduced to the form given by equation (92). When this structure, through relations (88) 



is fed into the evolution equation (85) one finds that Pgf^_^ — gf^_^P'^ = 0. This implies 
that the functions /^^ (^O)-^) are constants so that the solution gf{ko,kz) describes a free 
particle propagation without any quantum coherence. 

The above example may look trivial, but the important issue to note is that a precise 
information, or an ideal measurement of energy and momentum variables reduce the evo- 



lution equation (85) to a trivial description of a free particle propagation. So how does 
the usual density matrix picture with nontrivial quantum coherence emerge? The answer 
is that the coherence is possible only when the extrenous information about the state of 
the system is sufficiently inaccurate. Suppose now that our prior knowledge for example 
on the energy, momentum and spin variables can be described by some weight function 
]A^{ko,kz,s II k'Q,k'^,s'; z), where the primed variables are free and those without primes 
denote their known mean values. For example W could be a normal distribution in ko and 
kz, with variances ctq and az'- 

where we still took spin to be fixed and N is some normalization factor. (Other quantum 
numbers could of course be treated equally.) We now postulate that a physically observable 
density matrix can be defined in terms of the singular gf and the experiment related weight 

function as follows: 

■sr^ f dk' dk' 

pMko,kz,s;z) = y2 / ^^Wiko,kz,s\\k',,ki,s';z)g<iki,k',;z). (102) 

^ — ^ / Zn Zn 

s' 

"^This description serves the purpose for the present argument, and it would be the approach to be taken 
in many experimental situations. In reality, the role of the measurement is taken by the collision term that 
couples the system to the surroundings. This is actually the approach we shall take in the later publications 
where we will extend the present analysis to the case with collisions ^IM [T7] . 
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First note that our first example is easily described in this language, where we implicitly 
used a weight function which imposes strict ideal measurements of energy, momentum and 
spin of the state. As a second example, let us now assume that we have a complete ignorance 
on the direction of the momentum of the state, while we do have a precise information of 
the spin and the energy. Assuming that ko = u > this setting corresponds to a weight 
function 

W = 27i6{u-k'^)6s,s'- (103) 

It is now easy to see that the corresponding smeared out density matrix 



py^{uj,k,,s;z) = ^ 



dk' dk' 

° ' 2Tr5{uj-k'Q)5s,s'9f'ik'^,k'o;z) = p^i^; z) 



2n 27r 
obeys the standard evolution equation: 

isdzPs = P Ps - Ps P^ 



(104) 



(105) 



where 



P 



m* 



-m 

-OJ 



(106) 



is just the operator given in Eq. (83) with fco set to the externally imposed value fco = ^■ 



This equation is exact, and the singular structure plaguing the Eq. (85) has now been 



removed by integration so that Eq. (105) indeed is just the "normal" density matrix equation. 



capable of carrying information about coherent evolution. We stress again that this nontrivial 
structure emerged as a result of convoluting the (here rather the lack of the) available external 
information about the system on the definition of the physical density matrix. 



Because of the singular form of gf, the integration in (104) is trivial and the components 



of ps can be directly related to the on-shell functions /^ appearing in Eqs. (92) and (94): 

1 



pIl 



PRR 



Plr 



PRL 



-(- 



+ s)r- + 



s)fl 



I, UJ 



m 

TTl 

m* 



(/! + n 



9^k 



if-+f+) + 



(/r + ^/2 ) 



+ 2.;^i 

.rriRrs 

+ 2u;^i 



2 J' 

mi_ 

2uj^' 



(107) 



where /^ refer to functions /* {uj;z) with Sk^ = ±1 in Eq. (90) and /fg to fl2{uJ',z) in 



Eq. (94). Note that all components of p* mix terms from the mass- and coherence shells. 



One can extract the information about the particle numbers and coherence from ps at any 



point of the calculation by inverting the equations (107). Note that the four degrees of 



freedom encompassed by the combined mass-shell and coherence shells matches that of the 
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most general Herinitian 2x2 density matrix. Without coherence shells ps would contain only 
two degrees of freedom which is insufficient to describe any nontrivial quantum mixing. 
We conclude this subsection with comments related to the choice of the weight functions. 



First, it should be kept in mind that relations (107) between pij and fa are specific to the 



particular weight function Eq. (103). In principle, the weight connection could be something 
completely different, possibly encoding much more complicated structures of extrenous in- 
formation about the system. This information could be spatially dependent (as is the case 



in the example in section 5.2), or involve some partial, yet incomplete, information about a 
given quantum state. In any case, for any weight function there would always exist an in 
principle calculable relation connecting the two sets of variables. 



Finally, we point out that a weight function similar to ( 103 ) is actually the appropriate 
one to use for example for the interactions with a mass wall, encountered in the Electroweak 
baryogenesis problem. The reasoning is that in such case the only information one has 
about the system comes in the form of a set of conserved quantum numbers; in this case 
the energy, the momentum along the wall and the spin perpendicular to the wall. The 
density matrix can always be taken to be diagonal in conserved quantum numbers, but we 
can impose no extrenous constraint on a variable like k^ for example. Instead, one needs to 
introduce an explicit collision term which will give rise to damping terms that destroy the 
coherence generated by the interaction with the wall. This is of course the ultimate goal of 
our formalism, the results of which will be presented elsewhere |16j . 

To illustrate the use of the physical density matrices, and the corresponding choices of 
the appropriate weight functions, we next use our formalism to solve two known reflection 
problems. These examples will also further underline the neccessity of retaining the coherence 
solutions in the density matrix. 

5.2 Klein problem 



As our flrst example of the use of evolution equations (105), we shall consider a fermion 
reflecting off a step potential. This is of course a well known Klein problem, whose solution 
is known in the Dirac equation approach. For this problem we need to use the extended 



version of our master equation (25), with Ssing(M) = 7 V(m) where V{u) is the potential 



appearing in the usual Dirac equation. In addition we now take the mass to be a real 



constant so that equation (25) becomes: 

(^^^~m-7oF(u))G'<(M,i;)=0. (108) 

In the mixed representation one readily flnds 

{^+\p,-m--i^V{x)e-^^^^-^^)G<{k,x) = Q. (109) 
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Figure 4: Reflection from a step-like potential. Long arrows describe tlie directions of 
momenta of the in- and outgoing particles, and the short vector their direction of spin. 



Given this equation, we shall proceed with the analysis as in section (H). In the case of 
a step-potential the spatial gradients acting on V vanish everywhere except exactly at the 
potential wall, and within the wall the potential can be absorbed to the energy k^. Apart 
form the singular wall front the solution must then be of the form (95), where ko ^ k^ — V 



within the wall region. Moreover, interaction with a wall conserves the spin in z-direction. 
Thus, in the region I, shown in figure |4| the density matrix is a quantum mixture of incoming 
and outgoing (say) positive spin states with energy ko = u and in region II it describes a 
single outgoing s = 1 state with an effective energy ko = u — V. This information can be 
expressed in terms of a single 2;-dependent weight fuction W as followsjj 



W = 277 [9{z)6{ko -ij) + 9{-z)6{ko -u + V)] (5,,i. 



:ilO) 



The explicit form for the density matrix can now be derived from Eq. (102) 



p{z) = e{z)pi + e{^z)pii 



:nv. 



where pi^j can be directly (leaving out an overall half in all matrices) read off from Eq. ( 107). 
We find 



Pi = /i 



+ 1 



m 
k 



m 
k 



+ fl 



m 
k 



m 
k 

+ 1 



+ fl 



1 

m 



f, 



:ii2) 



where k = \Jufl — vrv^. For p\\ one replaces uj 



UJ 



V and k ^ q = ^^{uj - Vy 



m^ 



As 



usual there are two distinct possibilities depending on whether q is real or not. For a real q 
we expect a transmitted wave in the region II, but no incoming wave from the left, so that 
asymptotically f^ = 0. The function /^(-z) corresponds to the flux of reflected states and 



""In reality energy is conserved, but with this trick we get to account for the potential in a single weight 
function. 
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f-{z) to the flux of transmitted states, wliile functions f{2 (-2) give tlie colierence. Using 



Eqs. (107) we can write Eq. (105) dircectly in terms of the /-functions: 

sdji = 
sdzfi = -2kof2 

sdj2 = ^/i, (113) 

where k^ = to in the region I and ko = cu — V region II, and k^ = \/k^—m^. We see that the 
functions f^ are constants, and so we can normalize the incoming flux to one in the region 
I: /I = 1 and we find that f]^ = throughout. Moreover, the functions fi'2 are oscillatory 
for a real momentum. Since our boundary condition excludes asymptotic right moving states 
in the region II, also coherence functions must vanish there: fl\ = 0. In region I coherence 
is possible and one finds 

fl^{z) = Ai,2 cos{2kz) + 5i,2 sm{2kz) . (114) 



The coefficients A2 and B2 are related to Ai and Bi through equations (113) and the re- 
maining coefficients /^, f^, Ai and Bi are set by matching pi and pu at 2; = 0. One finds 
that 

f+ = ^ — , where x = — — ^ , (115) 

the flux is conserved: 

/+ + /" = 1 (116) 

and flnally the coherence solution is: 

/i(^) = ^:^T^cos(2M and fl{z) = -—dJliz) . (117) 

fc^g 1 + X 2uJ 

This is just the familiar result known from a Dirac equation approach [18], where f\. and 
f^ are identified with the usual refiection and transmission constants. In particular when 
l^ ^ we get q ^ k and x — > 1, so that f\_ — > 0, f^ -^ 1 and fl2 —* a.s expected. 

In case when q is imaginary, we cannot have any mass-shell solutions in the region II, 
so that /j? = 0. However, we can allow coherence solutions to be nonzero there, as long as 



they become asymptotically zero as z -^ —00. It is evident from Eq. (113) that when q is 
imaginary, appropriate exponentially decaying coherence solutions do exist. After a short 
calculation one finds the result with a complete refiection: /:j_ = 1 and with 

I = ^("^ ~/) cos{2kz) + ^ sm{2kz) (118) 

mV mV 

and f\{z) = —dzfl{z)/2uj. In the region II one has f]^ = and the coherence functions are: 

fl\z) = ^^^-^e^l^l^ (119) 

mV 
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Figure 5: Reflection from a potential wall due to spatially varying complex mass function. 
Note that in this example the mass is not related to energy in the same way as the potential 
V in the previous figure. 

and fl{z) = —dzfl{z)/2{u — V). The pure coherence in this case can be interpreted as 
describing a virtual pair of left moving state and right moving antistate (an anti-left mover). 
The lesson to be learned from this excercise is the necessity of including the k^ = 0-shell 
solutions in the mixture of the states; should we have dropped them, there would have been 
no consistent solution to the problem at all. This is not surprising, beause leaving out /i^2 
would physically correspond to making precise measurements of the momentum content of 
the state at all positions, arbitrarily close to the wall. Such measurements would significantly 
disturb and alter the quantum system under study. 

5.3 Reflection from a CP- violating mass wall 

As another reflection problem, we shall use our density matrix formalism to re-derive the 
CP-violating chiral flux from a wall created by a spatially varying complex mass function. 
This is the simplest example of a reflection problem relevant for electroweak baryogenesis, 
and it has been studied in the Dirac equation approach for example in references [HI [15] . 
The setup for the problem is depicted in figure [5] In this case the mass function is assumed 
to arise from the Higgs-mechanism, such that 



771(2) = y(f){z) 



(120) 



where y is a Yukawa coupling and (f){z) is some complex scalar field [15]. Deep in the unbroken 
phase the field vanishes, (j){z) ^ as 2; ^ 00, and all particles are massless. Far in the broken 



phase on the other hand, particles have a finite mass y(j){— 00) 



r7i_ 



whose phase can be 



chosen real. In the vicinity of the phase transition wall mass function is complex and spatially 
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varying. To be specific, we have used the following parametrization for 



l</>WI 



1 



ta.nh.{z / iw)) 



1 



arg( 



A^(l + tanh(z/0) 



:i2r 



where £^ is the width of the wall, and A^ is the total magnitude of the change of the phase of 



). This problem can be described by equations (83 86). The density matrix appropriate for 



the problem is set by the weight function (103), because all we know is that energy and spin 



are conserved quantities. Thus the equation of motion for our density matrix in chiral basis 



is just Eq. (105) where the density matrix components are given by Eq. (107). In practice 



it is more convenient to employ the Bloch-representation for the physical density matrix: 



Ps ^ T;i{9o) + {9") ■ ^) 



(g^) 



2tx 



9aikz,UJ;z). 



:i22l 



In this representation the equation of motion (105) becomes 



sd.ig'o) = 


~2mR{gl)~2mj{gl) 


SdM) = 


-2ko{g',)-2mM) 


SdM) = 


2ko{gl)-2mR{g^,) 


SdM) = 


0, 



:i23) 



and the variables (g^) are related to on-shell functions /, 



(a 



±, 1 or 2) as follows: 



Km. UJ U 



rriR 



I) 



fs) 



, (/- + /;) + /! 

-f^(/! + /;) + /2^ 
s{fi - /;) . 



(124) 



It is easy to solve equations (123) with a simple shooting algorithm. We take the initial 



conditions to be such that the incoming particle flux is normalized to unity, and that no flux 
is coming from the left. In terms of the on-shell functions /^ these conditions correspond 
toS 

^Indeed, note that with the normahzation (92), we find that the fermionic current is related to /|.: 

dfcod^fc. 



0-') = i^r^) 



E 



d-^fc 

(2^ 

dfcod^fc 



Tr[a''^G<] 
3 (/+ 



(27r)3 



s 



(2^) 



■(33 



(125) 



That is, fj. are to be interpreted as flux densities per unit energy and perpendicular momentum volume in 
the phase space. This is the appropriate interpretation for a problem where fc^ is not conserved globally. 
However, since on mass-shell and in the mean field limit d/co — f^zdfc^, the functions /|.'s can also be 
understood as local particle number densities per local unit 3-momentum. 
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Figure 6: a) Shown are the mass-shell functions f± corresponding to the left and right moving 
fluxes and functions /i 2, which encode the quantum coherence. We have taken s = 1 and 
g/|?7i_oo| = 0.088. b) The chiral density matrix components for the same solution. The wall 
parameterization is the same as in figure [7j 



/- 

/+ = /l = /2 



1, 

0, 



z 
z 



00 



-00 . 



(126) 



where we have set s = 1. In figure |6^ we plot the flux-functions f± along with the coherence 
functions /i^2 for the case where the asymptotic momentum to mass ratio is q/\Tn_ao\ = 
0.088. From Fig. (|6k) one can see that to the right from the wall the system is a coherent 
superposition of left and right moving states with opposite momenta, and the kz=^-s]ie\\ 
functions are oscillating coherently. In the broken phase however, all but the /.-function die 
off and the state soon becomes a pure transmitted left moving state. Note that the physical 
fiux-normalization condition /_|_(oo) + /_(— 00) = /_(oo) = 1 is satisfied by this solution. In 
Fig. rt6p) we plot the components of the chiral density matrix (107) for the same solution. 
The imaginary part of the plr goes to zero when 2 -^ —00, as a result of our choice that 
m becomes asymptotically real in the broken phase. The fact that the asymptotic state is 
an eigenstate of the effective Hamiltonian in the broken phase is seen in that p tends to 
a constant. Diagonal components of p become large in the broken phase. This is easy to 
understand, because Tr[p\ = (g^) is normalized to represent particle density (per unit energy 
and /C||-volume in the phase space, see footnote 6). However, since flux is conserved, a small 
local velocity Vz = kz/oJ enhances density in the broken phase by a factor of 1/vz- 

In figure [7] we show the particle-antiparticle fiux-asymmetry Aj_|_ = (/+ — f+)z=oo as a 
function of g/|m_oo| where q = y^co'^ — |?7i_ooP. In order to compute the density matrix 



for antiparticles we simply need to make the transformation m -^ m*. The characteristic 
peaked shape of the fiux-asymmetry arises as follows: for uj < |m_oo| particles cannot enter 
the broken phase at all, resulting in a complete reflection both in particle and in antiparticle 
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Figure 7: Shown is the current asymmetry of reflected states as a function of the asymptotic 
momentum to the mass ratio in the broken phase. For the wall width and the total change 
of the phase we used iy^ = 2 and A9 = — 1. 



sectors. When q is positive, the reflection amplitudes and hence the reflected asymmetry 
start to grow. This growth is cut for large q when the particle compton wavelength becomes 
shorter than the wall width and the particles start to behave classically, and the reflection 
amplitudes and hence also the asymmetry start to decline again (exponentially). These 
results were derived earlier for example ref. [T5] using the Dirac equation approach. We 
considered this case here partly to show how our formalism works in a nontrivial situation; 
in particular one sees again the necessity of including the k^ = 0-solutions into the definition 
of the full density matrix. Second, this application is precisely the one we wish to solve in 
completeness, including the collisions, in the EWBG problem 



5.4 Homogenous case, particle number in the early universe 

We introduced the concept of the weighted density matrix in the case of static planar sym- 
metric problem. The notion is of course more general and equally well adaptable to the 
homogenous problems. One interesting application of our formalism concerns the defini- 
tion of the particle number in the early universe. This problem was recently considered by 
Garbrecht et.al. in [3]. We shall now see how these results follow straightforwardly in the 
present formalism. Extension of this analysis to the case with interactions will be published 
elsewhere [T7| . 

In the homogenous case the spatial momenta k and the helicity are good quantum num- 
bers. Since the system is also isotropic, it is natural to consider density matrices of the 
form: 



pw{ko,\k\,h;t) = ^ 



dk'Q d^k' 



W{ko, \k\,h II k',, \k'l h' ■ t) g<{k'„ \k'\;t) . (127) 



27r (27r)3 
The physical situation considered in ref. |3j is that of a coherent particle production in 
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the early universe by an oscillating homogenous scalar field coupled to fermions. Unlike 
momentum and helicity, energy is not conserved and, having no prior information on it, we 
have to leave energy completely unspecified in the definition of the physical density matrix. 
Then the appropriate weight function is: 



Wi 



(27r) 



Anlk 



/|2 



Si\k\-\k'\)5h,h' 



(128) 



where the factor in front of deltas is introduced to normalize the 3-dimensional phase space 
density of states. It is now easy to see that the corresponding weighted density matrix 



Pmiko,\k\,h;t) 



E 



dk'o d^k' (2 



IT 



27r (27r)3 47r|fc'|2 



6i\k\-\k'\)6,,,,>g<ik',,\k'\;t) 



dfcp 

27r 



g<iko,\k\;t)= PhiHt) 



obeys the evolution equation: 



dtPh = -i[H,ph] 



(129) 



(130) 



which now is perfectly sensible, nonsingular equation where H is given by Eq. (39). If one 
introduces the Bloch-representation for the weighted density matrix: 



1 



Ph 



i{9o) + {f)-^), 



:i3i) 



one can write the equation (130) for the integrated components {g'^) = J ^ g'^ as 



dt{g1) 
dt{92) 
dtigl) 



2h\k\{g'^)~2mj{g',) 
-2h\k\{g1) - 2mn{g^) 
2mR{g!^) + 2mj{g1), 



[132) 



This equation is equivalent to Eqs. (31-40) in ref. [3|^ The authors of ref. |3] did not con- 
sider the constraint equations in their work, and averages {g'^) were introduced as unspecified 
moment functions, whose connection to the particle number had to be worked out using op- 
erator formalism and Bogolybov transformations. In our treatment the mass-shell functions 
/^ are directly related to the desired particle number densities. Indeed, for the fermionic 
current density we get: 



(^7» 



mjAo<i^i:Ji^M^i:l 



d^k 



M^ + f')- 



(133) 



^Note that the signs of our functions ((?2 3) differ from the corresponding ones in [5] because we define 
our 2-point functions differently (see footnote 1). 



31 



In the last step we used the exphcit form of (qq) in terms of /^'s; the complete set of 
expressions for all components is: 



(^o) = 


- /+ + /' 




(^l) = 






{92) - 


UJ 




{93) - 


- -h^-{f1-f1)+h(^f^- 

^ ^ k\ 


mi rh 



(134) 



According to the Feynman-Stueckelberg interpretation the actual phase-space particle num- 

ij:^ = f^{\k\) for fermions and n^^ 



ber densities are nr, = fl{\k\) for fermions and rir, = 1 — //^"^(Ifcl) for antifermions. Thus, 



using the inverse relations of Eqs. (134) we get for a given 3-momentum k and helicity h: 



i/,ir*i/i,» /^» , h^ \ i 



''^'^ - 2co 



h\k\{g',)+mn{g1)-mj{g^,)) + -{g',) (135) 



By setting a constraint Tr[p'*] = (g^) = 1 this reduces to the expression used to define the 
particle number in [3], apart from some sign conventions (see footnote 7). Similarly, the 
antiparticle number is found to be 

n^h = ^kh - i9o) + 1 • (136) 



kh 



Setting (g^) = 1 thus corresponds to assuming zero chemical potential: ((^q) = 1 ^ n 
n^^. In Fig. Mwe plot the particle number n^^ = fh+{\k\) and a function fc = {fli + fh2)^^'^) 
which measures the overall coherence between particles and antiparticles, in the case of a 
time dependent mass term, corresponding to an oscillating inflaton field during inflatonary 
preheating, introduced in ref. p]. We note that the generation of the particle number is 
highly coherent phenomenon with the amplitude of quantum coherence increasing with each 
oscillation period of the inflaton field. In ref. ^Jl\ we generalize our present formalism to the 
case with interactions and show how the interactions change the particle number production 
and how they introduce the quantum decoherence leading to eventual statistical ensemble 
of particles. 

6 Generalization to flavour mixing 

In this section we generalize our results to a case with several different flavours, i.e. when 
the mass function is replaced by an A^ x A^-matrix M. First note that this generalization 



has nothing to do with the chiral decompositions we made in arriving equations (30 31) in 



the time-dependent case or to equations (78 79 ) in the case with planar symmetry. The only 
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Figure 8: Shown is the fermion number density (thick hne) n^^ as a function of time r = 
uj^t for negative hehcity /i = — 1 in an inflatory preheating model defined in ref. [3]. For 
antifermions n 



kh 



Tij:^. Thin dotted line shows the overall amount of coherence between 



fermions and antifermions defined as /c = \/ f^i + fh2- Effects of inflaton oscillation are 
modelled by a varying mass term m{t) = (10 + 15 cos{2u!ff,t) + i sin(2c<j^t))|A;| where u^ = \k\ 
is the frequency of the oscillation. 



difference is that the mass operators in the generalized Hamiltonian (32) and momentum 



operators (80) are to be replaced with matrix operators. 

tdtg< = Hg< - g<H^ , 2k,g< = Hg< + g<H^ 
with 



H 



-h\k\ M 
ift h\k\ 



in the homogenous case and 



(137) 



(138) 



isd^gf = Pgl 



spj 



-2sk,gf = Pg; 



ifP^ 



(139) 



with 



P 



ko -M 
Mt -ko 



(140) 



in the planar symmetric case. Here we define mass operators M = M exp (— 19^ ■ dk) and 
M"'' = M't'exp (— 19*^ ■ dk), where ^^'■'^-derivatives always operate on the mass matrices, and 
9fc-derivatives on matrix functions gfi^. In the mean field limit the Hamiltonian operator H 
becomes a Hermitian local matrix operator H and the right hand sides of the kinetic and 



constraint equations in (137) become a commutator [H,g^] and an anticommutator {H,g^}, 
respectively. (Note that {MY 7^ M^ except in the mean field limit.) Equations (137 138) 



and (139 140) are completely general in the coUisionless limit. However, their interpretation 
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is complicated due to the same issues related to the singular spectral shell solutions to the 
constraint equations we found in the case of a scalar mass function. In what follows, we shall 
consider the effects of flavour mixing in the planar symmetric case. 

6.1 Planar symmetric case with flavour mixing 

As before, we shall constrain our analysis to the mean field limit M -^ M. Moreover, for 
the simplicity of notation, we shall take M to be Hermitian M = M'^ . Generalization to a 
non-restricted complex mass matrix is straightforward, but the non-Hermitian structure is 
not relevant for the qualitative issues we wish to discuss here. Introducing again a Bloch 



representation in chirality, we can write the constraint equations in (139) in a component 
form 



kogl + sk,g'Q-^-[M,g'^] 







0, (141) 

where M and g^, are NxN matrices in the flavour space. Let us again first study the case 



sk,g{ + ]^{M,gl} 



IS 



2k, 



[M, 



where k^ ^ 0. From the three last equations in (141) one obtains 

92 

9l 

g{ = -i-{M,^a + ^{M,[M,[M,^^]]} 



skz g s 

"k^^'^'Akjio 



2^^^'^°> + 8Pfco 



[M,[M,g'o]] 
1 



[U2) 



Putting these solutions back to the first equation in ( 141 ) then gives the spectral equation: 

1 



ikt-ki-M')g^, + -[M 



2 „si 



16A;2 



[M^[M^^^]] =0. 



(143) 



In order to carry the analysis further, we need to go to the basis where the mass matrix is 
diagonal. Since M is assumed to be Hermitian, there is a unitary matrix U such that 

tud = UMU^ (144) 

is diagonal. Correspondingly, we denote the density matrix in the diagonal basis by: 

9ha = Ug^U^ . (145) 
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The first observation to be made is that the two last terms in the equation ( 143 ) are purely 
off-diagonal in the mass eigenbasis: 









(146) 



where Am?- 
form: 



m'jjj — m'l)-. Using the results (146), we can write (143) in the component 



fVQ 



K - Mf, - 



16/c2 



JoOJij 



0. 



(147) 



where M^- 

'■J 



(ml + m|) 



)/2r| Let us first consider the diagonal entries in the equation (147). 
Because /S.m\ = it is immediately clear that diagonal equations have solutions analogous 
to the solutions found in the case with a scalar mass function. The dispersion relation is 



fczii — 31 \ / fCg 



m 



Di 



and the corresponding spectral solution for the density matrix elements: 



{9 



DO) a 



"^^fls. ih, z) -— 6{k, - Sk, Jkl - m 



"Di) 



'1481 



(149) 



Moreover, it is easy to see that for the diagonal elements the equations (142) reduce to 



equations (88) and similarly that the diagonal parts of the evolution equations reduce to the 



equivalent expressions in the scalar mass case. We then obtain the following solutions for 
the chiral structure of the diagonal matrix elements: 



{9 



sDJii 



1 



(9 



DO) a 



1 skz/koii 
mDi/koii 



moi/koii 
1 + skjkoii 



(150) 



which, as before, is seen to describe free propagation of a given helicity mass eigenstate. 
Before moving on to discuss the off-diagonal constraint equations on mass-shells, let us now 
find out if any kz = 0-solutions might be left out by our previous derivation. Setting kz = 0, 
the constraints become 



'2ko9Do-{'mD,9Di} = 

'2ko9m-^['rnD,9D2\ = 

['mD,9m] = 

{mD,9m} = 0' 



The last of these equations immediately implies that g'^{k 
constraint on the third line on the other hand only sets g^j^^. 



#j 



(151) 

0) = 0. The commutator 
: at A;^ = 0, but leaves the 



^This dispersion relation was found in ref. [5D], but the physical content of the small /c^-branch as one 
corresponding to the quantum coherence was not realized by the the authors of that paper. 
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diagonal elements of gf^Q arbitrary. Taking commutators and anticommutators of first and 
second equations, respectively, with respect to niD similarly leads to constraints (^Iji^. , = 0, 
but again leaves g^j^^ 2 arbitrary. Restricting now to diagonal part of the first constraint on 



(151), we easily find 

rriD 



9 m,. = ^9m,.^ (1^2) 

dX kz = Q. That is, we find that all off-diagonal components of g^)^ must vanish on shell 
kz = 0. Moreover, also the diagonal elements of gf^^ vanish while the components g£,i2 



remain arbitrary and the nonzero diagonal elements of (/f^g are given by (152). From our 
previous results it is obvious that these solutions encode the information of the quantum 
coherence between the mass eigenstates of opposite helicity and momentum. The different 
on-shell solutions with their interpretations for the mass eigenstates are shown in the figure 
M for the case of 2 x 2-fiavour mixing. 



Let us now turn to the off-diagonal dispersion relations in Eq. (147). These solutions de- 
scribe the quantum coherence between different mass eigenstates. The most striking feature 
about the off-diagonal dispersion relation is that for large ko it has two distinct solutions; 
one for k^ close to the diagonal shell momentum and another one close to k^ = 0. Moreover, 
no solutions for the dispersion relation exist for ko < max(?7ii5j,?7i£ij). The interpretation of 



these solutions is easier when we rewrite the dispersion relation (147) in a different form: 



ihh = ±1 [^jkl-ml^ ± ^kl - ml)i . (153) 

In this form we now explicitly indicate also that the momentum shell depends on the off- 



diagonal entry in question. The first signs in Eq. (153) refer just to the direction of the mo- 
mentum, while the second two signs refer to the high- and low-/c^ branches of each continuous 
dispersion curve, to the right and to the left from the minimum set by fco = max(r7i£)j, niDj) 
respectively. This structure is depicted in Fig. [9] in the 2x2-mixing case. The off-diagonal 
momentum shells thus correspond to the mean momenta of the mixing mass eigenstates. We 
can now make the following physical interpretations: 



• The large momentum solutions (with plus sign inside the parenthesis in Eq.(153) carry 
the information of the quantum coherence related to flavour mixing between different 
mass eigenstates of equal helicities moving in the same direction. 

This is of course just the quantum coherence phenomenon that is relevant for example for 
neutrino-oscillations. A triplet of such potentially mixing states with positive kz is shown 
by a circle in Fig. [9j Second, 



The small momentum solutions (with minus sign inside the parenthesis in Eq.(153) 
carry the information of the quantum coherence related to flavour mixing between 
different mass eigenstates of opposite helicities moving in opposite directions. 
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Figure 9: The complete set of dispersion relations corresponding to Eqs. (147) or (153) 



These solutions are the analog of the k^ = 0-solutions found to host the quantum coherence 
relevant for the reflection problem in the case of a scalar mass and for the diagonal entries 
in the case of a mass matrix. At this point one might appreciate the fact that we did 
not obtain any kz = 0-shell contributions for off-diagonal terms above; these solutions were 



already included in the "on-shell" dispersion relations (147). In fact, Eq. (153) actually 
contains also the kz = 0-solution for the diagonal entries, as this is what the small-fc^ 
solution collapses to when the masses men and moj are equal. This solution is simply 
hidden in the form (147) of the dispersion relation when one first takes the limit Amf, -^ 0. 
In Fig. |9] we show the complete set of dispersion relations for this problem. The thick 
dash-dotted red and green parabolas show the diagonal mass-shell dispersion relations and 
thin dash-dotted lines at A;^ = show the coherence shells corresponding to the diagonal 
solutions. Triplets of (red and green) dots on the dash-dotted lines show particular sets of 
phase space elements involved in the potential diagonal mixing. The thick solid blue line 
shows the off-diagonal dispersion relation and the triplets of (brown) dots show shells that 
are involved in the off-diagonal flavour mixing between states moving to opposite directions. 
Finally, the set of three black dots enclosed in a circle show a triplet of shells involved in 
the usual flavour mixing between states moving to a same direction. We have used different 
values of ko in different cases just for the sake of clarity of presentation. In reality, in a 



case with no information on kz (like in the reflection cases we considered in sections 5.2 5.3) 
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there are nine different Icinematical shells with sixteen unknown functions, contributing to 
a physical density matrix that describes the most general mixing phenomena. This is the 
setting one expects to encounter (but with a generic, non-Hermitian mass matrix) in the 
case of a chargino reflection off a phase transition wall in an application of our formalism to 
the Electroweak baryogenesis. 

6.1.1 Evolution equation for a neutrino beam 

Having now explained the physical significance of different shells, we move on to briefly 
discuss the dynamical equations in the case of a free propagation. In the case of reflection 
problem the new aspect is the need to account for the mixing, as induced by locally varying 
unitary operator U{z). We shall not consider this problem in its full generality in this paper. 
Instead, let us now ignore the small fc^-branches entirely and concentrate on flavour mixing 
between states of (nearly) equal momentum. That is, we are assuming that we have a good 
resolution of the direction and the magnitude of the momentum, but not good enough to 
separate the different flavour mixing shells from each other. For simplicity we will assume 
perfect information of /cq here. This is of course just the case of interest in case of a neutrino 
beam, or a neutrino flux moving in a given direction, such as reactor or solar neutrinos. 
Now, the matrix evolution equations written in the component form for the physcial density 
matrix are: 

sdM)-mA9l)] = 

sdM) + {M,{92)} = 
sdM)-'2ko{gl) + {M,{g',)} = 

sdM) + 2h{g2)+i[M,{93)] = 0. (154) 



Using the integrated form of the third equation in (141) we can eliminate ((7|) from the 
equation for {g^), which becomes simply 

^5.(^^) = [M^(|-)]^^[M^(^^)], (155) 

where k corresponds to the mean momentum as given by the observational accuracy for- 
mally encoded in a weight function W for the problem. It is in fact easy to verify that all 



components {g^) obey an identical evolution equation, such that (155) actually describes the 



flavour evolution of a freely propagating state of a given deflnite spin and momentum, but 



with explicit flavour mixing. Of course the spin, which no more shows explicitly in Eq.(155) 
can be replaced by helicity, by associating the 2;-axis with the direction of motion of the 
particle. 



Let us make a couple of remarks on the solution (155). First, it is just the analog of the 



free particle motion found in sections 4-5; either for a scalar mass function, or more generally 
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for all diagonal elements the commutator in Eq. (155) vanishes. Second, we can interpret 



the equation (155) as arising from a Hamiltonian form: 



id,p=[H,p], (156) 

if one makes the formal identification H = \^k^ + M^ ~ fc + M"^ /2k. In this case equations 



(155) and (156) are identical to first order in M'^/k'^. Note however that the equation (155) 
is in fact exact (apart from the averaging in the last step) and its definition exactly encodes 
the amount of the information (momentum resolution) on the system. It should be obvious, 
that at the limit when the resolution encoded by W becomes accurate enough to separate 
different shells in the flavour mixing triplets shown in Fig. |9} the density matrix becomes 



trivial and the commutator term vanishes in Eq. (155), and the solution is reduced to a free 
propagation of a fixed mass, momentum and helicity eigenstate. Note that our formalism 
allows one to consider also the intermediate cases where one has partial information on the 
flavour content of the state, such that this information affects, but does not stop completely 
the mixing and oscillation pattern. 

7 Interacting fields 

The goal of this paper was to set up the density matrix formalism for treating quantum 
coherence phenomena in classical backgrounds. Eventually we wish to extend these methods 
to include cases with collisions. We will not pursue this goal further here, apart from a 
qualitative discussion of the elements needed in the derivation. For the time dependent 
case the generalization is actually quite straightforward, and we will present our complete 
results in a companion paper [17]. Most results found in this paper were derived from 
the free coUisionless equation (0) for the correlation function iGc{u,v) = (Tq [^(u)-?/'(f )]), 
defined on a complex Keldysh time-path shown in Fig. [T] More generally, in the presence of 
interactions, Gc{x,y) obeys the contour Schwinger-Dyson equation: 

Gciu, v) = G^u, v)+ I dSi I d^Z2 Gl{u, ^i)Sc(;2i, Z2)Gc{z2, v) , (157) 

where S^ is a self-energy functional, which in general depends on higher order Green's 
functions of the theory. This dependence eventually leads a hierarchy of coupled equations 
involving all possible Green's functions. The practical usefulness of the Schwinger-Dyson 
formalism arises from the fact that in many applications one can truncate this hierarchy to 
the lowest order by some reasonable approximation to Sc which only involves the 2-point 
functions. In the weak coupling limit for example, it is natural to do this by substituting 
all higher than 2-point functions by their perturbative expressions. A recent review that 
discusses the evaluation of Sc can be found in reference [H]. 
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Figure 10: A generic form of a Schwinger-Keldysh equation for the 2-point function G* 



Equation (157) is formally expressed in Fig. 10, where the thin lines correspond to the 



free particle (tree level) propagator G^, and the thick lines to the full propagator Gq, and 



the filled ellipsis represents the self-energy function Sc. Multiplying Eq. (157) by the inverse 



of the free particle propagator (6*°) ^ and integrating over the connecting variable zi one 
finds 



{i P^ ~m*PL - mPR)Gc{u, v) 



(5c(m° - v'')5%u -v)+ dSSc(M, z)Gciz, v), (158) 



where Sc{u^ 



is a contour time delta-function. Here we assumed the free Lagrangian 



of the form Eq. ([2| . It thus appears formally obvious how the formalism can be extended 
to the case with collisions; one merely needs to evaluate the appropriate function Ec and 
proceed in the derivation as described in this paper. There are several obstacles on the way 
to a set of equations that can be solved in practice however. The crucial issue turns out to 
be finding an approximate way to treat the phase space of the interacting system in a way 
that retains the notion of a single particle excitations. This can be done in a meaningful 
way in the so-called quasiparticle and the mean field limits (or up to first order in gradients 



for fermions). Taking these limits, the mixed representation equation (19) becomes just 



(^ + 2 ^^' ~^^ ~ irh^T 



j:R)G'^{x,k) =Ccoii{x,k) 



(159) 



where Hr is the real part of the (retarded) self energy function and Ccoii is the collision 
integral. In the thermal equilibrium approximation Ccow can always be written as 



a 



coll 



-^T (G< - G<) 



(160) 



where G^ is given by Eq. (62) and F is the usual thermal collision rate. The quasiparticle 



approximation is familiar from thermal field theory [2T]. In the present context it corre- 
sponds to neglecting all terms arising from Ccow in the constraint equations. Under these 



assumptions, equations (159) can be shown to support a spectral solution for the phase space 



with (quasiparticle) mass and coherence shells, similar to the ones described in this paper. 
Given this structure to the phase space, one can define physical density matrices as weighted 
integrals, and compute how the collisions affect the particle distribution functions related 
to the various mass and coherence shells. The resulting formalism can be used to describe 
for example the effects of collisions on coherent particle production in the early universe 
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and approach to thermal equihbrium including quantitative account of the emergence of 
decoherence [TTj. 

Let us finally note that deriving an interacting theory for the static problems, is some- 
what more subtle, since the usual CTP-formalism leads to Green's functions that describe 
correlations that vanish at temporal infinities. This is consistent with the usual definition 
of the asymptotically free vacuum states for the theory using temporal infinity. This is not 
the appropriate limit for the static reflection problems, where one rather would like to see 
correlations vanish at spatial infinities. Correspondingly one would like to define the vacuum 
states of the theory at spatial infinities, and develop a scattering formalism relating vacua 
and states at different spatial rather than temporal infinities [16j. 

8 Conclusions and outlook 

In this work we have derived quantum kinetic transport equations for fermionic systems 
including non-local quantum coherence. A crucial observation leading to our formalism was 
the finding that in cases where the full translational invariance is lost, the free fermionic 
2-point correlation functions G^'^{k,X) have, in addition to the usual mass-shell solutions 
k"^ — |?7ip = 0, new solutions living on shells ko = (homogenous case) or kz = (planar 
symmetric static case). These solutions were identified as carrying the information about 
the quantum coherence between particle-antiparticle pairs of same helicities and opposite 
momenta in the former case, and between incoming and reflecting waves of equal spin in 
the latter case. Another crucial element was the definition of a physical density matrix as a 
convolution of the singular phase space density matrix with a weight function encoding the 
amount of extrenous information (or the quantitative measure of the lack of it) about the 
state of the system. We illustrated the use of the formalism with several examples including 
reflection problems and definition of a particle number in the early universe during coherent 
particle production. We considered also the case of several mixing fermion fields and showed 
how the usual evolution equation for fiavour mixing neutrino system arises from the singular 
phase space structures and an appropriate weight function. Finally we have outlined how our 
formalism can be extended to the case with interactions, with a slight technical reservation 
concerning the spatially varying problem. 

Let us finally comment on our choice to limit the discussion to the mean field limit. This 
seems somewhat contradictory, since quantum effects become more important when the rate 
of change in the background gets larger in comparision with the wave length or the frequency 
of the probe. Also, we got exactly the correct answers to our refiection calculations despite 
the mean field limit assumption. Understanding these apparent paradoxes begins from the 
observation that (at least for sufficiently smooth weight functions), the integrated evolution 
equations always have the same form as in the mean field limit, since all derivative corrections 
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to them are reduced to vanishing surface terms. So, the only thing that does get changed 
by gradients is the connetion between the averaged-out density matrix elements and the 
mass- and coherence shell distribution functions. Indeed, our particle numbers and fluxes 
are just mean field approximations to the full quantum system. However, these connections 
do become exact in the case of the Klein problem everywhere outside the potential step, 
and in the case of smooth wall at spatial infinity. This is of course why the formally mean 
field quantities in these cases give exact results for asymptotic currents. Similarly, in the 
coherent particle production case, our particle number and coherence functions provide an 
approximation to the full quantum phase space, that becomes exact when inflaton oscillation 
stops and mass becomes a constant. 

One might wonder if these considerations render our results to be only of academic 
interest. This is obviously not so: first the complicated structure of the phase space is not 
optional; it is there. One cannot just ignore the constraint equations and concentrate to 
the integrated form of the evolution equations. We have shown that in the mean field limit 
this structure is singular and allows a particularily transparent picture for separation to 
quantum coherence and mass-shell degrees of freedom. It is true that beyond mean field 
limit, the singular structure is lost. In the case of fermions this occurrs at the second order 
in gradients, while the first order can be computed within spectral limit and it gives rise 
to corrections that lead to the semiclassical effects discussed in refs. [H [3 |9]. Even then 
the mean field limit can provide a good approximation to the phase space, capturing the 
most important features of the quantum evolution. Second, from the practical point of view, 
the singular shell structure for mass- and coherence shells is crucial when the formalism is 
extended to include interactions. Indeed, the entire success of the current approach relies 
on ones ability to find a spectral approximation to the dynamical and kinematical phase 
space of the system; only then can we compute the collision terms explicitly and describe 
the evolution of the coherence and particle numbers on these shells in a tractable manner. 

As we have pointed out in many occasions, this paper is merely setting up the basic 
formalism which will be extended elsewhere to include decoherence P^ [T7] and then applied 
to various problems of interest in cosmology. This formalism will be crucial in particular 
to reliably compute the quantum reflection contribution to the baryon number production 
during the electroweak phase transition [191 ttOl US] . It will also be possible to use it to study 
the effect of collisions on the coherent particle production [17] . We believe that the formal- 
ism could, and also will provide to be useful in other applications beyond the immediate 
application to the cosmology. 
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